Awaking the vacuum in relativistic stars 
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Void of any inherent structure in classical physics, the vacuum has revealed to be incredibly 
crowded with all sorts of processes in relativistic quantum physics. Yet, its direct effects are usually 
so subtle that its structure remains almost as evasive as in classical physics. Here, in contrast, we 
report on the discovery of a novel effect according to which the vacuum is compelled to play an 
unexpected central role in an astrophysical context. We show that the formation of relativistic stars 
may lead the vacuum energy density of a quantum field to an exponential growth. The vacuum- 
driven evolution which would then follow may lead to unexpected implications for astrophysics, 
while the observation of stable neutron-star configurations may teach us much on the field content 
of our Universe. 
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Born from the efforts to formulate a consistent rela- 
tivistic version of quantum mechanics, the successful for- 
malism of quantum field theory (QFT) has unveiled an 
incredibly rich structure for what was thought in clas- 
sical physics to be the most uninteresting of the states: 
the vacuum. From the early Dirac sea of negative-energy 
states to the picture of virtual particles constantly be- 
ing created and annihilated, the vacuum has acquired 
conceptual importance for a consistent description of na- 
ture. However, from the observational point of view, its 
existence continues to be almost as evasive as in classi- 
cal physics, demanding carefully designed experiments in 
order to detect its subtle effects. (The Casimir effect 
according to which two (or more) objects experience a 
force between them solely due to the fact that their pres- 
ence changes the vacuum energy in a position-dependent 
way, is such an example, successfully observed in labo- 
ratory. See, e.g., Refs. @, Q.) Particularly interesting is 
the fact that the vacuum, being a dynamical entity, gravi- 
tates, which using the concepts of general relativity (GR) 
means that it affects and is affected by the geometry of 
the spacetime. This rests at the root of the Hawking ef- 
fect [j, [B| , according to which black holes should emit a 
thermal bath of particles. Although this semiclassical- 
gravity effect of "black hole evaporation" has become a 
rare landmark for those who seek a complete quantum 
gravity theory, it is virtually unobservable for realistic as- 
trophysical black holes. Here, in contrast, we investigate 
possible astrophysical implications of a recently-proposed 
scmiclassical-gravity mechanism [6] by which the vacuum 
of free quantum fields is forced to play a central role. We 
focus on the context of relativistic stars and show that 
the formation of (realistic) compact objects may disturb 
the vacuum of a quantum field in a way which leads its 
energy density to an exponential growth. Eventually, this 
vacuum energy should take control over the evolution of 
the background spacetime, determining the ultimate fate 



of the relativistic star and possibly leading to unexpected 
implications for astrophysics. On the other hand, the 
observation of stable neutron-star configurations may be 
used to rule out the existence of fields for which the effect 
should have been triggered. 

Vacuum energy in gravitational fields: The task 
of calculating "vacuum energies" of quantum fields in 
gravitational backgrounds is one which is properly ad- 
dressed (though not easily performed in general) only by 
using techniques of quantum field theory in curved space- 
times (QFTCS), which is the minimal and most conser- 
vative way of merging the successful formalisms of QFT 
and GR |7H12| . Born more than four decades ago 0, HI , 
QFTCS has proven its value by providing the proper con- 
text in which legitimate "low-energy" (i.e., below Planck 
scale) quantum gravitational phenomena, like the black 
hole evaporation, could be unveiled. Therefore, before 
presenting the arguments which lead to our main result, 
let us give the reader a brief overview on how to obtain 
the vacuum energy of a quantum field in a curved back- 
ground using the canonical approach to QFTCS. In a 
nutshell, it goes as follows: (i) one solves the generally co- 
variant field equation for a complete set of normal modes 
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} of the field $ in a given fixed background 



geometry with appropriate boundary/initial conditions 
which express the choice of the (Fock) Hilbert space of 
interest; (ii) one (formally) expands the field operator <J> 
using these normal modes with the operator coefficients 
of the positive- norm (negative-norm U& ) modes 

being interpreted as the annihilation (creation) opera- 
tors for that particular mode; (iii) then one (formally) 
constructs the stress-energy-tensor operator T^ v out of 
its classical expression by substituting $ by and calcu- 
lates its expectation value in the state of interest (usually 
the vacuum of the chosen Fock Hilbert space). Since 
is a quadratic expression in <£> (which should be viewed 



as an operator- valued distribution), some regularization 
and renormalization procedure must be applied in step 
(iii) in order to give rise to an expectation value which is 
free from ultraviolet divergences. 

Because the background spacetime has to be fixed al- 
ready in step (i), the procedure outlined above is ex- 
pected to give meaningful results whenever the backre- 
action of the expectation value of T^ v on the background 
spacetime (through the semiclassical Einstein equations) 
can be neglected. Notwithstanding, failure to comply 
with this condition can be used to expose situations 
where semiclassical-gravity effects become crucial. This 
is the strategy we follow here. 

Disturbing the vacuum during the formation of 
a star: For the sake of simplicity, we follow Ref. @ 
and consider a real scalar field $ with mass m > 
and arbitrary coupling £ with the Ricci scalar curva- 
ture R. The field equation is the usual Klein-Gordon 
equation with the additional scalar-curvature coupling: 
(-□ + m 2 + £i?) $ = 0. (We adopt natural units in 
which h = c = 1 unless stated otherwise.) We consider 
the case in which matter initially scattered throughout 
space with very low density eventually collapses to form a 
static (and stable) spherically symmetric compact object 
(according to classical GR). Therefore, the background 
metric can be taken to be asymptotically flat in the past 
and asymptotically static and spherically symmetric in 
the future: 
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where / = f(x) > and r — r(x) > are functions of 
the coordinate \ alone, with f(x) — > 1 and r(x)/x ~~ * 1 
for x oo, and dr/dx > so that no trapped light-like 
surface is present. 

The natural in- vacuum state is the one associated with 

(u~ )* which in the 



the normal modes and tti ' 
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asymptotic past take the form of the usual flat-space sta- 



tionary modes: 
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where k € M 3 and := + m 2 is interpreted as the 
energy of the mode. Since the spacetime is also asymp- 
totically static in the future, after the compact object 
is formed, there is a natural out-vacuum state, which is 
the one associated with a different set of normal modes 
v usiix an d v< wili = ( v miu)* wm ch in the asymptotic future 
assume the stationary form 



(+) future 
zolfj. 
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, / ) are the usual 



where (/ = 0, 1, 2, ... and n = —I, . 
spherical-harmonic functions and F m i are solutions of 
Schrodinger-type equations 



(-d 2 /d x 2 



V 



(I) 



F m i = vo 2 F„ 
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satisfying usual boundary conditions and normalization. 
Here, w is the energy of the mode (which no longer needs 
to be greater than m) according to static observers in the 
asymptotic future which are far away from the compact 
object. The effective potential v}g = V}g (x) is given by 

V$ =f[m 2 +ZR + 1(1 + l)/r 2 ] + (l/r)d 2 r/d X 2 . (5) 

The functions f(x) and r(x) can be related to the mat- 
ter/energy distribution of the compact object through 
the Einstein equations of GR. For instance, assuming the 
compact object to be made of perfect fluid (assumption 
which will be used throughout this article for simplicity) , 
the effective potential can be put into the form 
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where p — p(x) is the matter/energy density of the com- 
pact object and p = p(x) represents its averaged value 
up to x- P{X) ■= 3M(x)/(47rr(x) 3 ), with M( X ) being 
the mass of the object up to x- (Recall also that, ac- 
cording to GR, R — 8irG(p — 3p), where p is the pressure 
which sustains the configuration of the compact object 
and G is Newton's gravitational constant.) 

In general the in- and out- vacuum do not coincide, 
which means that particles are created due to the change 
in the gravitational background. However, this parti- 
cle creation is usually negligible. In contrast, as recently 
pointed out in a quite general context [6| , there is another 
possible influence of the gravitational background on the 
quantum vacuum: well-behaved gravitational fields may 
trigger an exponential increase of the vacuum energy den- 
sity leading to a vacuum-dominated scenario. Applied to 
the context of compact objects, the gravitational fields 
which are entitled to "awake" the vacuum energy are 
the ones for which the effective potential in Eq. © gets 
sufficiently negative somewhere in the asymptotic future 
so that Eq. (f4]) possesses well-behaved solutions F^i for 
vo 2 = — Q 2 < 0. In this case, additional, non-stationary, 

"ml and w ail = ( w nii)* with the asymp- 



totic behavior 
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are necessary in order to expand the field in the asymp- 
totic future. This implies that at least some of the 
in-modes eventually go through a phase of exponential 
growth, which leads to an exponential increase of the vac- 
uum expectation value of <t 2 . This is then reflected in an 
exponential growth of the vacuum expectation value of 
f^v (see Eqs. (10-13) of Ref @). (It is important to note 
that this asymptotic divergence has nothing to do with 
high-frequency modes, so that renormalization plays no 
important role in this analysis.) As a consequence, no 
matter what the value of the classical energy density of 



2 



the compact object is, eventually the vacuum energy den- 
sity (and pressure, and momentum density) of the field 
will dominate and take control over the spacetime evolu- 
tion. 

Sharp transition to vacuum dominance in rela- 
tivistic stars: It is interesting to note that in principle 
the vacuum-dominance effect could be triggered for mas- 
sive fields with m 2 <C R ~ Gp (so that a region where 
IQj is negative might exist for values of |£| not too large). 
However, recovering the units in which these quantities 
are usually expressed, this condition implies 



m 2 /(3.5 x 1Q- 12 eV) 2 
p/(10 14 g/cm 3 ) 



< 1. 



(8) 



Therefore, considering that relativistic stars can reach 
densities p ~ 10 14-17 g/cm 3 [l3j |. the most "natural" 
mass satisfying this is m = 0, which we shall assume in 
what follows. 




FIG. 1: Diagram showing the values of the ratio M/R s of 
uniform-density compact objects which trigger the exponen- 
tial growth of the vacuum energy density of massless scalar 
fields with coupling £ (dark-gray region). The vertical dashed 
line indicates the conformal-coupling value £ = 1/6. The 
black region represents the values of M/R 3 for which no equi- 
librium can be classically obtained. 

It may come as a surprise to verify that the simplest 
distribution possible for an idealized compact object, the 
uniform energy-density distribution, is already entitled 
to awake the vacuum energy of quantum fields. As- 
suming p = p = po, the Tolman-Oppenheimer-Volkoff 
(TOV) equation of relativistic hydrostatic equilibrium 
can be analytically integrated to obtain the pressure dis- 
tribution p and then the function / (see, e.g., Ref. 14|). 



Then, substituting these into Eq. (J5J) , a numerical search 
for "bound" solutions of Eq. flU with uj 2 < (and 
1 = 0, which is the most promising case) can be per- 
formed for several values of the radius R s and the mass 
M = 4irpoRs/3 of the compact object. In this uniform- 
density case, the existence of such "bound" solutions de- 
pends only on the ratio M/R s and the value of £. The 
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FIG. 2: Same as diagram in Fig. [T] but now for compact ob- 
jects with parabolic density profile. The light-gray region rep- 
resents values of M/R s for which equilibrium configurations 
are classically unstable for most realistic equations of state 
(see, e.g., Ref. [IB]). 



results of this search are shown in Fig. [TJ We see that a 
massless field with any value of£>l/6or£<— 2 can 
have its vacuum energy density exponentially amplified 
for some range of M/R s , 

This idealized uniform-density case serves to illustrate 
that well-behaved background geometries can indeed in- 
duce the vacuum-dominance effect for fields with appro- 
priate masses and scalar-curvature couplings. This, how- 
ever, is not enough if we want to explore possible ob- 
servational implications. For that matter, more realistic 
density and pressure profiles for the compact object have 
to be used. Although there is a certain amount of un- 
certainty on the correct equation of state for matter at 
densities as high as 10 15-17 g/cm 3 , the density profile 
of realistic neutron stars with masses 1.5M© ~ 2.2M 
can be closely approximated by a parabolic radial depen- 
dence, p = p c (l - r 2 /R 2 s ) 0, where p c = 15M/(8nR%) 
is the central density. In this case, the TOV equation 
can also be analytically integrated to obtain the pres- 
sure distribution p and the background metric [16| , and 
a numerical search for "bound" solutions of Eq. (fj| with 
vo 2 < (and I — 0) can be performed, as in the previous 
case. The existence of "bound" solutions again depends 
only on the ratio M/R s and the value of £. The results of 
this search are shown in Fig. [2] Note that Figs. [T]and[2] 
are very similar but it is worthwhile to point out that in 
this more realistic case the effect is triggered for smaller 
values of M/R s . In particular, values of M/R s for which 
most realistic equations of state are causal [15| , and which 
represents classically stable configurations, can awake the 
vacuum of fields with proper couplings. 

The time scale which governs the exponential growth 
of the vacuum energy density, energy flow and pressure 
is given by (which usually is of order |V^g \~ 1 ^ 2 )- 
In fact, the dominant contribution to the vacuum energy 
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density in the asymptotic future can be estimated from 
Ref. 0: 

future /"» 7 /_\ 2Qt i t-»3 

p v ~ flh(r)e zut I R 6 S 

~ fe(f)exp ^ /(1 °" 5s) V°' 62g/Cm3 (9) 
Air; exp (^ /(10 km) J R4 j {1Q km)4 , W 

where h(f) is a dimensionless function of f := r/i? s and 
we have used that for a bound solution f2 ~ 1/R S - The 
exact values assumed by h(f) depend on the details of 
the spacetime evolution between the two static asymp- 
totic regions, but for f ~ 1 its order is not far from unit 
(though it is exponentially damped for large values of f). 
Therefore, for a compact object with radius R s 10 km 
whose effective potential happens to possess a bound so- 
lution, it would take only a few milliseconds (according 
to far away observers) for the vacuum energy density to 
become overwhelmingly dominant over classical matter 
densities as high as 10 14 ~ 17 g/cm 3 . This suggests that 
the transition from classical to vacuum dominance can 
be very sharp. In spite of this, we emphasize that the to- 
tal vacuum energy is conserved (as it should in this static 
scenario where backreaction is not considered) due to the 
fact that the volume integral of h{f) vanishes. 

Final Comments: Although we have been talking 
about "vacuum dominance" of quantum fields, we would 
like to clarify that the appearance of exponentially grow- 
ing modes already occurs in the classical level. This 
is evident from our argument that the modes given in 
Eq. (J7J are necessary in order to expand an arbitrary 
(classical) field configuration. Therefore, the very same 
situations which trigger the exponential growth of the 
vacuum energy density will most likely lead to an ex- 
ponential growth of the classical energy density stored 
in an initially non-vanishing field configuration. Here, 
however, we have opted to focus on the quantum version 
of the effect since it is triggered even in the absence of 
any initial classical perturbation, which makes its con- 
sequences more robust than the ones coming from its 
classical counterpart. The existence of initial classical 
perturbations can change the details but not the conclu- 
sion about whether or not the compact object will face 
the exponential instability, which is the main issue here. 

The vacuum-awakening effect unveils a quite interest- 
ing interplay between semiclassical gravity and observa- 
tional astrophysics. Concerning field theorists, the obser- 
vation of some relativistic stars may be used to rule out 
the existence of certain fields in Nature. For instance, 
according to our Fig. O the observation of a stable and 
approximately spinless cold neutron star with mass M 
and radius R s (with an approximately parabolic density 
profile) rules out the existence of massless scalar fields 
with couplings £ for which (£,M/R 3 ) lie in the dark- 
gray region. Considering that 95% of the energy content 
of the Universe is unknown, finding ways of testing the 
existence of free fields (for which there would be no rea- 



son to assume an initial classical perturbation) is very 
welcome. Finally, concerning astrophysicists, the awak- 
ening of the vacuum energy of certain fields may deter- 
mine the ultimate fate of some relativistic stars. This 
will eventually depend on how the vacuum energy back- 
reacts on the spacetime, as ruled by the semiclassical 
Einstein equations. The vacuum-driven evolution may 
lead, e.g., to (i) a complete collapse, hiding the negative- 
effective-potential region inside an event horizon, or (ii) 
an explosion, ejecting mass and bringing the effective po- 
tential back to a profile admitting no bound solution (or 
possibly a combination of both). Only a detailed nu- 
merical analysis, beyond the scope of this article, can 
decide what will happen once backreaction is taken into 
account. Either way, we would have an unusual situation 
where the vacuum would be responsible for an event of 
astrophysical proportions. 
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